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The true- and eccentric-anomaly parametrizations of the Kepler motion are gen- 
eralized to quasiperiodic orbits by considering perturbations of the radial part of 
kinetic energy as a series in the negative powers of the orbital radius. A toolbox of 
methods for averaging observables in terms of the energy E and angular momen- 
tum L is developed. A broad range of systems governed by the generic Brumberg 
force, as well as recent applications of the theory of gravitational radiation in- 
volve integrals over a period of motion. These integrals are evaluated by using the 
residue theorem. It is shown that the pole of the integrand is located in the origin 
and that under certain circumstances an additional pole emerges. 

In our previous works we have analyzed the motion of a compact binary 
system with rotating components in order to provide realistic signal templates for 
the intcrfcromctric gravitational-wave observatories. We have calculated the cumu- 
lative radiation reaction effects on the motion of the spinning binary. In obtaining 
secular changes one has to evaluate integrals of type 



where T is the radial period, n an integer and u> a function with no manifest radial 
dependence. For example T is determined by n — —2, while averaging instanta- 
neous radiative losses of quantities characterizing the orbit of a binary system one 
calculates this type of integrals with n = 1, . . . , 7. During the computations it is 
convenient to use the generalized true- and eccentric-anomaly parameters defined 
as follows. 

The radial equation of a perturbed Kepler problem is considered as 



where the small coefficients ipi couple the perturbations to the system. The conser- 
vation of energy E and magnitude of the orbital angular momentum L are assumed 
to hold from symmetry considerations. A wide class of perturbed Kepler motions 
accomplish these conditions, for instance systems described by the Brumberg^ force 
and the dynamics of spinning binaries. 

The generalized true- and eccentric-anomaly parametrizations are introduced 
by the relations 
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where 7 and re are positive constants. The scale of the parameters is adjusted 
to the turning points r(0) = r min , r(w) = r max , which are solutions of the radial 
equation f = 0. To first order one has p — 2 unphysical roots (even if real, they are 
much smaller than the radius r) and r m ; n and r max , which differ from the Keplerian 
turning points in linear terms in ipi. With the above properties the familiar form 
of the parametrizations is obtained : 

2 1 + cos y 1 — cos v „ 1 + cos £ 1 — cos £ ... 
- = + , 2r = + 3^ . (4) 

T Tmin '"max r min T'max 

The most remarkable properties of these parametrizations are summarized in 
the following theorems: 

Theorem 1: The integral J r~^ 2+T ^dt , n > is given by the residue located 
in the origin of the plane of the complex true-anomaly parameter z — e lx . 

Theorem 2: The integral J r~( 2+n ^dt with n < is given by the sum of the 
residues in the origin of the complex parameter plane w = e 1 ^ and [for perturbative 
terms with powers i > 2 — n of r in the radial equation (TJPJ at 

1 /2 

/ mu 2 - J-2riEL 2 \ 

\mfi 2 + ^-2n,EL 2 ) 

The above statements are proved by passing to the true- and eccentric-anomaly 
parameters in the integral and examining this expression in linear order of the 
perturbations (pi. In doing so one can specify the possible to functions which will 
not destroy the advantageous property of the integrand. It turns out, that u can 
be any polynomial in sine and cosine of the parameters. In particular, it can have 
the form 

U3 = ojq + <p(i>,r) , (5) 

where lu is constant and the perturbing term ip(tp, f) is a polynomial in sin^, cos tp 
and r. To first order the expressions of the azimuthal angle ip and r in terms of the 
true- and eccentric-anomaly parameters are needed to Keplerian order. 

The two types of generalized Keplerian parametrizations studied here are ca- 
pable of describing a wide variety of systems ([2]). They do not give a complete 
parametrization of the motion, but are suitable for handling radial effects. When 
using them and the residue theorem the evaluation of various integrals ([T]) turns 
out to be a simple task. They have be en t ested in computing the secular radia- 
tive changes of a compact binary systenP^. In our view the solution of the r adial 
motion presented here is at least a viable candidate amongst other approaches 
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